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Abstract. The baryon production mechanism at the intermediate pr (2 - 5 GeV/c) at RHIC is still not 
well understood. The beam energy scan data in Cu+Cu and Au+Au systems at RHIC may provide us 
a further insight on the origin of the baryon anomaly and its evolution as a function of t/snn- In 2005 
RHIC physics program, the PHENIX experiment accumulated the first intensive low beam energy data in 
Cu+Cu collisions. We present the preliminary results of identified charged hadron spectra in Cu+Cu at 
ytSNN = 22.5 and 62.4 GeV using the PHENIX detector. The centrality and beam energy dependences of 
(anti)proton to pion ratios and the nuclear modification factors for charged pions and (anti)protons are 
presented. 

PACS. 25.75.Dw 



1 Introduction 

One of the most surprising observations at RHIC is a par- 
ticle type dependence of hadron yield suppression at the 
intermediate transverse momentum pt (2-5 GeV/c) 012]- 
In Au+Au central collisions at ^/snn — 200 GeV, yields 
for mesons are largely suppressed j^j with respect to the 
yields in proton-proton collisions in the intermediate Pt, 
while those for baryons are not suppressed, and show a 
binary nucleon-nucleon collision (N co u) scaling. 

In d+Au experiment at RHIC, there is also a parti- 
s' 1 cle type dependence in the nuclear modification factor 
' Raa EH- The Raa for protons in d+Au is larger than that 
for pions and kaons. This observations can be understood 
by Cronin effect It is found, however, that the par- 

ticle species dependence in d+Au is not large enough to 
account for the absence of suppression for protons as seen 
in Au+Au central collisions at RHIC. This phenomena in 
Au+Au central collisions is called "Baryon Anomaly at 
RHIC" . To explain the baryon anomaly, many theoretical 
models have been proposed, such as quark recombination 
models [7], and a hybrid model of hydrodynamics with a 
color glass condensate (CGC) and jet quenching All of 
these models are able to reproduce the experimental data 
qualitatively for both pions and protons. 

On the elliptic flow (112) measurements, there is a also 
clear particle species dependence [0]. The number of con- 
stituent quark scaling of v-i shows a clear universal curve 
regardless of the particle species, and these scaling prop- 
erties support the quark recombination picture for the 
hadron production at RHIC. 

In order to test the quark recombination picture, <f> me- 
son is one of the ideal test particles, because the mass is 
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similar to protons', but it's a meson particle. If the hydro- 
dynamical collective flow is a dominant source of baryon 
anomaly, (j) meson's Raa should follow the data for pro- 
tons (no suppression) due to its heavy mass. Thanks to 
the high statistics data in Au+Au collisions at ^/smn = 
200 GeV taken in 2004 RHIC run, it is confirmed that 
the Raa for cj> mesons behaves like a meson |10|. and its 
vi follows the same universal quark number scaling curve 
of all other particle species Giving the fact that </> 
mesons behave like a " meson" , it is now widely believed 
that the quark recombination is one of the main mecha- 
nisms for the hadronization process at the intermediate 
Pt in central Au+Au collisions at RHIC. 

Now the key questions are; where the onset of the 
baryon anomaly is, and how they evolve as a function of 
beam energy. In order to answer these questions, the lower 
energy beam data in Cu+Cu collisions were taken during 
the 2005 RHIC run by the PHENIX experiment. In this 
paper, we present the preliminary results of single particle 
Pt spectra for identified charged particles in Cu+Cu colli- 
sions at ^/snn = 22.5 GeV and 62.4 GeV measured by the 
PHENIX experiment. The beam energy and centrality de- 
pendences of the particle ratios and nuclear modification 
factors are presented and discussed. 



2 PID charged particle spectra analysis 

The low beam energy data sets in Cu+Cu have been taken 
during the 2005 RHIC run. We analyzed 5.2 M minimum 
bias events for 22.5 GeV Cu+Cu, and 66.5 M minimum 
bias events for 62.4 GeV Cu+Cu data set. We constructed 
the tracks within a collision vertex ± 30 cm from the cen- 
ter of the spectrometer. 
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Cu+Cu -> p+X @ \Jsj^ = 22.5 GeV 
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Fig. 1. Centrality dependence of pr spectra for protons in Fig. 2. Centrality dependence of pr spectra for protons in 
Cu+Cu collisions at ^/sjvjv = 22.5 GeV. No feeddown correc- Cu+Cu collisions at ^/sjvjv = 62.4 GeV. No feeddown correc- 



tion is applied. 



tion is applied. 



Charged particle tracks are reconstructed by the Drift 
Chamber (DC) based on a combinatorial Hough trans- 
form. The first layer of Pad Chamber (PCI) is used to 
measure the position of the hit in the longitudinal direc- 
tion (along the beam axis). When combined with the loca- 
tion of the collision vertex along the beam axis (from the 
Beam-Beam Counter(BBC)), the PCI hit gives the po- 
lar angle of the track. Only tracks with valid information 
from both the DC and PCI are used in the analysis. In 
order to associate a track with a hit on the high resolution 
Time-of- Flight detector (TOF), the track is projected to 
its expected hit location on the TOF. Tracks are required 
to have a hit on the TOF within ±2cr of the expected hit 
location in both the azimuthal and beam directions. The 
charged particle identification (PID) is performed by us- 
ing the combination of three measurements: time-of-flight 
from the BBC and the TOF, momentum from the DC, 
and flight path length from the collision vertex point to 
the hit position on the TOF. 

The geometrical acceptance and the in-flight decay are 
corrected by the GEANT based Monte Carlo simulation. 
No occupancy correction (detector inefficiency correction 
due to the high track densities) is applied. The expected 
occupancy correction is below 4% in Cu+Cu collisions 
62.4 GeV at the most central events and less for periph- 
eral collisions, the correction is negligible. No weak decay 
feeddown correction is applied for proton and antiproton 
spectra in this analysis. 



The centrality determination in Cu+Cu collisions 22.5 
GeV is purely based on the PCI hit multiplicity in or- 
der to avoid an effect of the spectator nucleons. Since at 
y/SNN — 22 GeV the beam rapidity gap is quite narrower 
than those at the higher beam energy, the BBC (located 
in 3.0 < |?7| < 3.9) is able to see the spectator nucleons. 
The PCI multiplicity distribution with the BBC z ver- 
tex correction is subdivided into four centrality bins. For 
Cu+Cu 62.4 GeV data, we use the charge distribution 
measured by the BBC. We obtain the number of binary 
collisions (N co u) and the number of participant nucleons 
{Npart) for each data set and each centrality bin by the 
Glauber Monte Carlo calculation in PHENIX. 



3 Results 

We present the following preliminary results in Cu+Cu 
collisions at y / ijvjv = 22.5 GeV and 62.4 GeV, measured 
by the PHENIX experiment. (1) the identified charged 
particle pr spectra, (2) particle ratios (p/n + , p/n~) as a 
function of px for each centrality class, and (3) the nuclear 
modification factor Raa for pions and (anti)protons as a 
function of px ■ 

3.1 pt spectra 

Fig. n and Fig. [21 show the centrality dependences of pr 
spectra for protons in Cu+Cu collisions at ^/sJjn — 22.5 
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Fig. 3. Centrality dependence of pr spectra for p in Cu+Cu Fig. 4. Centrality dependence of pr spectra for p in Cu+Cu 
collisions at ^/snn = 22.5 GeV. No feeddown correction is collisions at ^/sjvjv = 62.4 GeV. No feeddown correction is 



applied. 



applied. 



GeV and 62.4 GeV, respectively. The pt spectra for an- 
tiprotons are shown in Fig. [3] and No weak decay feed- 
down correction is applied for both protons and antipro- 
tons spectra and both beam energies. We measured the 
spectra in four centrality bins (0-10%, 10-30%, 30-60%, 
60-88 or 60-100%) and in the minimum bias events. The 
Pt spectra for each centrality bin for 7r and K are also 
measured at both beam energies (see H2J). 

3.2 Particle ratios (p/n + and p/ir~) 

Fig. shows the centrality dependence of p/ir + and p/tt~ 
ratios as a function of pt in Cu+Cu collisions at ^/snn = 
22.5 GeV. The similar plots for Cu+Cu at 62.4 GeV are 
shown in Fig. [B] 

In 22.5 GeV Cu+Cu data, p/ir + ratio is significantly 
larger than those in the higher beam energy in Au+Au pQ, 
and it is increasing as a function of collision centrality. 
This large p/7r + ratios can be understood by the influ- 
ence of the incoming nucleons from the beams. For p/ir~ 
ratios, (almost) no centrality dependence is observed, and 
the ratios are constant at the value of 0.3 - 0.4 at pt = 2.0 
GeV/c, which is close to the expected values in the frag- 
mentation and the data in p+p collisions. 

In 62.4 GeV Cu+Cu data, the influence of incoming 
protons is smaller than the one in Cu+Cu collisions at 
22.5 GeV, as seen in p/tt + ratios in Fig. The centrality 
dependence of p/tt~ is seen. In the most central collisions, 
p/tt~ ratio is about 0.6 at pr = 2 GeV/c, and in the 



peripheral collisions, the value becomes consistent with 
that for p+p collisions. 

Fig. shows the beam energy dependence of p/ir + and 
p/tt~ ratios in Cu+Cu collisions from ^/snn = 22.5 GeV 
to 200 GeV . There is a clear beam energy dependence 
in Cu+Cu from 22.5 GeV to 200 GeV, i.e. p/ir+ (p/tt") 
decreases (increases) as a function of ^snn, respectively. 

For the study of the baryon anomaly at the lower beam 
energies at RHIC, p/tt~ ratios would be a good measure, 
because antiprotons are "produced particles" , but mea- 
sured protons are the mixture of produced particles and 
the incoming protons from the beams. The absence of cen- 
trality dependence inp/7r _ ratio in Cu+Cu 22.5 GeV may 
suggest that there is no baryon anomaly at this beam en- 
ergy. To conclude this observation, the high statistics data 
with the heavier collisions system like Au+Au at around 
22.5 GeV is necessary in the future RHIC run, and the 
results should compare with the existing data in Pb+Pb 
collisions at the top SPS energy (^/snn — 17.3 GeV) at 
high p T HSj . 



3.3 Nuclear modification factor: Raa 

To obtain the nuclear modification factors at the lower 
energy, we used the pr spectra at mid-rapidity in p+p 
collisions at y/s — 23 GeV and 63 GeV for kaons and 
(anti)protons as the reference spectra, measured by the 
British-Scandinavian Collaboration jTHj. Those data are 
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Fig. 7. Beam energy dependence of p/ir + (left) and p/n (right) ratios as a function of pr in Cu+Cu collisions from ^/sjvjv 
22.5 GeV to 200 GeV. The data for 200 GeV Cu+Cu is taken from the PHENIX preliminary data [TiJ . 



p/7t* Cu+Cu 22.5 GeV 



PHENIX Preliminary 



p/jr Cu+Cu 22.5 GeV 



V" ' 

2 - PHENIX Preliminary 


1 i i i i 


i i i i l j. 


a- 


•o 


0-10% J 




■ .71 


10-30% 


5 j- 


AA 


30-60% ~z 


4 - 


▼ 


60-100% -: 


2- 
1- 

8 r 




Min. Bins_ 


5 r 


% 

i , , . , 


- I 



Pt [GeV/c] 



Pt [GeV/c] 



Fig. 5. Centrality dependence of p/7r + (left) and p/tt~ (right) 
ratios as a function of pr in Cu+Cu collisions at ^/snn = 22.5 
GeV. No feeddown correction is applied. 



pliC Cu+Cu 62.4 GeV 

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

2t- PHENIX Preliminary 



p/jr Cu+Cu 62.4 GeV 

"p 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

2- PHENIX Preliminary 
1.8L tO 0-10% 

10-30% 

18 f A A 30-60% " 

1.4 r TO 60-88% ■ 

i 91 Min. Bias. 



Pt [GeV/c] 



Pt [GeV/c] 



Fig. 6. Centrality dependence of p/tt + (left) and p/ir~ (right) 
ratios as a function of pr in Cu+Cu collisions at ^/snn = 62.4 
GeV. No feeddown correction is applied. 



fitted by the empirical functional form (see 22 m de- 
tail). For pions, we use the parameterization as described 
in 221- The nuclear modification factor Raa is denned as 
the following equation: 



Raa(pt) 



{l/N e A ^)d 2 N AA /dp T dy 
(N co u)/ai n p el x (Pa^/dpTdy' 



(1) 



is the average Glauber nuclear 



where the (N coll ) / a* p n p el 
overlap function, (TauAu)- 

The Raa for pions and (anti)protons are shown in 
Fig. El for 22.5 GeV Cu+Cu, and in Fig. H for 62.4 GeV 
Cu+Cu. The centrality selection for Raa is 0-10% central 
collisions for both beam energies. The values of the average 
number of collisions used here are (N co u) =140.7 (ai ne i = 
32.2 mb) for 22.5 GeV Cu+Cu, and (N coll ) =152.3 {a inel = 
35.6 mb) for 62.4 GeV Cu+Cu. The systematic error on 
Raa consists of; (1) < N co u > uncertainty (one a error 
as shown in dotted-lines above and below Raa = 1), (2) 
combined systematic errors from p+p reference data and 
Cu+Cu spectra (not shown in the figures). The statistical 
errors are shown in the error bars on each data point. 

The data show that there is no suppression on charged 
pions for both 22.5 and 62.4 GeV Cu+Cu central colli- 
sions. For proton's Raa, a clear enhancement is seen and 
can be attribute to the effect of the incoming beam nu- 
cleons. The value of proton's Raa at 62.4 GeV is slightly 
smaller than that at 22.5 GeV. Raa for antiprotons is al- 
most unity at the intermediate px for both 22.5 and 62.4 
GeV, and is very similar to the one for pions. 

In 2006, the PHENIX experiment has successfully col- 
lected the new reference p+p data at y / s7v"/v = 62.4 GeV. 
The precision of Raa measurements at 62.4 GeV for both 
Cu+Cu and Au+Au is expected to be improved by using 
this new p+p data set. 
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Fig. 8. Nuclear modification factor Raa (0-10%) for charged 
pions and (anti)protons in Cu+Cu collisions at ^/sjvjv = 22.5 
GeV. The positively charged particles are shown in the left 
panel, and the negatively charged particles are shown in the 
right panel. 
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Fig. 9. Nuclear modification factor Raa (0-10%) for charged 
pions and (anti)protons in Cu+Cu collisions at -y/Ijvjv = 62.4 
GeV. The positively charged particles are shown in the left 
panel, and the negatively charged particles are shown in the 
right panel. 



4 Conclusions 

In summary, we have measured pt spectra for tt^, , p, 
p in Cu+Cu collisions at ^snn — 22.5 and 62.4 GeV. 

In 22.5 GeV Cu+Cu data, we observed a larger p/ir + 
ratio compared to those at the higher beam energy in 
Au+Au, increasing as a function of collision centrality. 
For p/ir~ ratio, almost no centrality dependence is seen, 
and the ratio is about 0.3 - 0.4 at pt — 2.0 GeV/c, 
which is close to the value in p+p collisions. In 62.4 GeV 
Cu+Cu data, p/n + ratios are reduced, compared to those 
in Cu+Cu 22.5 GeV, and the centrality dependence of 
pl~K~ is seen. There is a clear beam energy dependence 
on those ratios in Cu+Cu from 22.5 GeV to 200 GeV, i.e. 
p/ir + (p/ir~) decreases (increases) as a function of y/SNN, 
respectively. The observed p/k~ ratio may suggest there 
is no baryon anomaly at 22.5 GeV Cu+Cu collisions. To 
make this point clearer, a further data set (a high statis- 
tics low beam energy data in heavy collision system) is 
necessary in the future RHIC run. 

For the Raa, no suppression is observed for charged 
pions for both 22.5 and 62.4 GeV in Cu+Cu central col- 
lisions. For protons Raa, there is a clear enhancement, 
which can be attributed to the incoming beam nucleons. 
The Raa in 62.4 GeV is slightly smaller than that in 22 
GeV. The antiprotons Raa is almost unity at the inter- 
mediate pt region for both beam energies. 
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